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We use a renormalized mean field theory to investigate the superconducting properties of un- 
derdoped cuprates embedded with overdoped or metaUic regions that carry excess dopants. The 
overdoped regions are considered, within two different models, first as stripes of mesoscopic size 
larger than the coherence length and then as point impurities. In the former case we compute 
the temperature dependent superfiuid stiffness by solving Bogoliubov de Gennes equations within 
the slave boson mean field theory. We average over stripes of different orientations to obtain an 
isotropic result. To compute the superfiuid stiffness in the model with point impurities we resort 
to a diagrammatic expansion in the impurity concentration (to first order) and their strength (up 
to second order). We find analytic expressions for the disorder averaged superfiuid stiffness and 
the critical temperature. For both types of inhomogeneity we find increased superfiuid stiffness, 
and for a wide range of doping enhancement of Tc relative to a homogeneously underdoped system. 
Remarkably, in the case of microscopic impurities we find that the maximal Tc can be significantly 
increased compared to Tc at optimal doping of a pure system. 
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I. INTRODUCTION 

Local probes of the cuprate superconductors reveal 
signatures of electronic inhomogeneity both at the mi- 
croscopic scales of lattice constants and at somewhat 
larger mesoscopic scales '^' . The inhomogeneity is gen- 
erally seen as separated regions with either a large or a 
small gap, which have been attributed to local variations 
in the doping level with respect to the half filled Mott 
insulator ' . 

Experiments which probe global properties indicate 
that the average doping level has two direct effects on the 
superconducting properties. First, the pairing gap is seen 
to decrease with hole doping away from half filling""'". 
Second, the superfiuid stiffness extracted from penetra- 
tion depth measurements, increases with doping . 
This interplay between two energy scales relevant to 
superconductivity is thought to give rise to the dome 
shaped dependence of on hole doping'-. Doping inho- 
mogeneity is therefore expected to lead to spatial modu- 
lations of the pairing amplitude along with variations of 
the charge carrier density. 

In this paper we shall investigate how inhomogeneity 
in the doping level affects global superconducting proper- 
ties of the material. Specifically we address the effect of 
inhomogeneity on the temperature dependent thermody- 
namic stiffness and, ultimately, on the transition temper- 
ature. To this end we employ a semi-phenomenological 
model of a d-wave superconductor that takes into account 
the the proximity to the Mott insulator through a strong 
on-site repulsion. Furthermore we consider various scales 
of inhomogeneities, ranging from the microscopic scale of 
a lattice constant to mesoscopic scales, somewhat larger 
than the coherence length (see Fig 1). An important 
question for practical applications is whether the transi- 
tion temperature can be enhanced significantly by judi- 



cious design of the inhomogeneity. The idea is to gain 
from an optimal combination of large pairing gap in the 
low doping regions and large carrier density in the highly 
doped ones^ '. 

Enhancement of Tc due to a similar mechanism was 
predicted in cuprate heterostructures composed of an un- 
derdoped superconducting layer coupled to an overdoped 
metallic one.^"'^'^'' The underdoped layer induces a prox- 
imity gap in the overdoped layer, which then contributes 
to the zero temperature phase stiffness of the system and 
considerably enhances it compared with the suppressed 
stiffness of the underdoped layer. On the other hand, 
the d-wave proximity gap which is induced on the metal- 
lic layer is small, and thus results in a sharp reduction 
of the stiffness with the temperature^''. We found in 
Ref. IG that the combined effect can in principle lead to 
enhancement of Tc compared with an optimally doped 
layer. However to attain such enhancement the cou- 
pling between layers needs to be much larger than the 
realistic coupling between the cooper-oxide planes. It is 
therefore unlikely that these simplified models provide 
a satisfactory explanation for the Tc enhancement ob- 
served in various experiments on heterostructures.^ '^^^ 
However, if there is doping inhomogeneity within a plane 
the coupling between the overdoped and underdoped re- 
gions would naturally be large since they are connected 
by the in-plane rather than the c-axis tunneling. As we 
shall see this situation can indeed give rise to enhance- 
ment of the maximal critical temperature compared to a 
pure system. 

Specific kinds of in-plane inhomogeneity and their ef- 
fect on superconductivity have been previously investi- 
gated theoretically. For example a weak-coupling BCS 
theory of the attractive Hubbard model showed that Tc 
can be enhanced by periodic modulations of the weak 
attraction."" A density matrix rcnormalization group 
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(DMRG) study of the repulsive Hubbard model on a two 
leg ladder showed that modulations of the hopping ma- 
trix element along the ladder can enhance the pairing 
correlations and thereby possibly increase the Tc of a cou- 
pled ladder system.^ ' A direct study of the two dimen- 
sional Hubbard model using contractor rcnormalization 
(CORE) also indicated that there is an optimal modula- 
tion of the hopping matrix clement, which maximizes the 
pairing correlations."" Finally, dynamical mean field and 
cluster Monte Carlo calculations find increased pairing 
gap, and possibly T^, in a state with charge modulation 
near 1/8 doping. ^'^'-^^ 

The above studies focus on the effect of periodic com- 
mensurate charge modulations on the pairing order pa- 
rameter. We complement and extend the analysis in sev- 
eral ways. First, we use an effective theory, amenable 
to analytic treatment that allows to identify the physical 
origin of the various effects. Second we compute the tem- 
perature dependent superfluid stiffness, which at least in 
the underdoped cuprates is a more complete measure of 
superconductivity than the pairing amplitude and allows 
us to directly estimate Tc. Third, in addition to the stripe 
model treated in previous work we also consider random 
doping variations, which appear to be the more generic 
situation in samples of doping above 1/8. Both for the 
stripe model and the random inhomogeneity we asses the 
possibility of enhancing Tc by tuning the magnitude of 
characteristic doping modulations and their length scale. 

We implement the inhomogeneity in the form of inclu- 
sions of a highly overdoped phase, already in the metal- 
lic regime, embedded in a background of underdoped or 
optimally doped material. The case of mcsoscopic in- 
homogeneity, where the metallic inclusions are of the 
size of the superconducting coherence length or larger 
is sketched in Fig. 1(a). This is treated within an ef- 
fective stripe model of the metallic regions, where we 
average over stripe orientations to obtain an isotropic 
macroscopic stiffness. Another case we consider, is where 
the metallic regions are much smaller than the coherence 
length and are modeled as point impurities. This case is 
depicted in Fig. 1(b). 

In both cases we include the crucial effects of strong 
coulomb repulsion and of the d-wave symmetry of the 
order parameter. The former is the reason for the low 
superfluid density ps at low doping, while the second 
is responsible for the linear suppression of ps with T 
at low temperatures." ' These effects are taken into ac- 
count within a slave boson mean field theory of the 
t — J model."'"'"' Furthermore, wc include Fcrmi-liquid- 
like corrections phenomenologically, to the description of 
low energy quasiparticles."''"'^" 

Regardless of the model for the metallic regions we 
find an increase of the zero temperature stiffness and for 
a wide range of doping levels, also higher critical temper- 
ature compared to the pure system with the same average 
doping. Furthermore, in the case of microscopic impu- 
rities we even predict that a higher Tc can be attained 
even compared to the maximal Tc at optimal doping of 




FIG. 1: (a) An illustration of a mesoscopic-scale inhomo- 
geneous layer. The typical size of the metallic regions is 
equal or larger to the superconducting coherence length, (b) 
Microscopic-scale inhomogeneous layer, the metallic regions 
are point-like impurities. 



the pure system. 

The paper is structured as follows: In Sec. II we give a 
general overview of the models used, of the assumptions 
that underlie our choice of models, and of the main re- 
sults obtained in the different regimes. Section III gives 
a detailed treatment of a model representing mcsoscopic 
inhomogeneity, while in section IV we consider a model 
with point impurities. Section V is a summary and dis- 
cussion of the results. 



II. OVERVIEW 

In this section we introduce the framework for treating 
the inhomogeneous cuprate layer within a slave boson 
mean field theory. We describe the essential ingredients 
of the theory for the case of mesoscopic inhomogeneity 
as well as for point impurities. Finally we summarize the 
main results that are derived in detail in later sections. 

In order to describe doping inhomogeneity in cuprate 
materials we make use of models that can account for 
the effects of doping of the Mott insulating parent com- 
pound. A simple theoretical framework that captures 
many of the important effects is the renormalized mean 
field theory (RMFT)-'' or slave boson mean field theory 
(SBMFT)-' of the t - J Hamiltonian, 

HtJ = ~Pg ^»J-4cj-^G + h.C + Y, J^J^^^3■ (1) 

ija (ij) 

Here = Atf^/U is the super-exchange interaction. 

Si = clgdss'Cis' and Pq — Hi(l — ni^riii) implements the 
Gutzwiller constraint, which prohibits double occupancy 
of sites. 

The standard mean field treatment of the t — J model 
includes two approximations. The first is to account for 
the projection only through renormalization of the hop- 
ping tij — >■ gijtij, while working in the full rather than 
the projected Hilbcrt space-''. The second approximation 
consists of a standard decoupling of the quartic term in 
both the Fock and BCS channels. The resulting mean 



3 



field Hamiltonian is given by 

HmF = (gljt + Xij)4aCja + h.C-Y^ tJ-iC^aCta 

i.j,(T i^s 

+ Y.^^A^]i + ^-'' (2) 

where Xi] ^ S-^ij E<r(4<TCj<T)/4, A^j = 3Jy(4^c]^ - 
c|j^c^^)/4 and gij are doping dependent renormalization 
factors that account for the effect of the no-double- 
occupancy constraint. In a uniform system of doping 
p, g — 2p/ (1 +p), Xij — X for nearest neighboring i, j, 
and l^i^i^x = —^i,i+y = A such that the pairing has a 
dr^2_y2 symmetry. 

The mean field theory of the t — J model captures the 
crucial fact that the zero temperature superfluid stiff- 
ness of undcrdoped cuprates scales linearly with the hole 
doping, po oc p.^"'-'' It also accounts for the d-wave sym- 
metry of the gap that gives rise to a low energy quasipar- 
ticlc spectrum of the form = [(v/k||)^ -I- {va^±)^]^^^- 
This form of the spectrum explains the observed linear 
reduction of the superfluid stiffness with temperature, 
Ps{T) = Pa- boT with bo = 2\og2{2y/2Ztf / {-KVfVA)-^'^ 
However, the mean field theory does not give the cor- 
rect value of Z . This can be viewed as a Fermi liquid 
correction that may be strongly renormalized at low en- 
ergies due to quasi-particle interactions not included in 
the mean field theory.-"""'" Therefore Z is best taken 
as a phenomenological parameter to be extracted from 
experiments.-"''^^ 

In this paper we extend the analysis of the stiffness and 
the critical temperature to the case of an inhomogeneous 
system. Specifically we describe an underdoped system 
in the bulk (0.1 < pi < 0.15) embedded with highly 
overdoped metallic regions. We consider two regimes of 
inhomogeneity as illustrated in Fig. 1. First is when the 
metallic regions are of the order or larger than the super- 
conducting coherence length and second when they are 
of the order of one lattice constant. As discussed in the 
introduction a pertinent question we wish to address is 
whether such inhomogeneity can lead to enhanced Tc- 

A. Mesoscopic Inhomogeneity 

In the first model, described in Sec. Ill, we assume a 
2D mixture of a superconducting underdoped phase (of 
doping pi) and an extremely overdoped metallic phase (of 
doping p2 > 0.3). The doping level varies considerably 
only across a length scale of the order of the coherence 
length ^ ^ Vf / A., which is typically around 5 lattice spac- 
ings, such that the 2D regions are of intermediate size > ^ 
as depicted in Fig. 1(a). This scenario is reminiscent of 
various experiments that find gap variations on a similar 
scale, of 5 — lOnm^""'"' . 

We model this system as a mixture of striped do- 
mains, each one with alternating underdoped and over- 
doped stripes along the x or y direction, such that on 



a macroscopic scale the system is fourfold rotationally 
invariant [see Fig. 2(a)]. This allows us to obtain an 
expression for the superfluid stiffness of the entire sys- 
tem. The superconducting stripes are described by the 
t — J Hamiltonian and the metallic stripes are modeled 
by free fermions. We vary the widths of the stripes in 
order to explore the superconducting properties in vari- 
ous geometries. To calculate the critical temperature of 
the inhomogeneous mixture, we solve self consistently the 
Bogoliubov de Gennes equations for Hamiltonian (2) al- 
lowing for position dependent gij , A^j and Xij ■ We derive 
a general formula for the superfluid stiffness Ps{T) of a 
striped superconductor in terms of response kernels that 
can be directly calculated from the Bogoliubov de Gennes 
solution [see Eqns. (7), (11)]. We then use the Kosterlitz- 
Thouless criterion Ps{Tc) = 2T^ln to determine of the 
mixed system. 

We show that there exist optimal configurations which 
allow for an enhanced zero temperature superfluid stiff- 
ness in the inhomogeneously doped layer, compared with 
the homogeneous superconducting one. This is a con- 
sequence of proximity effect that leads to a gap in the 
metallic regions. The metallic regions, having a large 
density of charge carriers, can then contribute signifi- 
cantly to the superffuid stiffness of the inhomogeneous 
layer at T = 0. On the other hand, since the proxim- 
ity gap is much smaller than the original superconduct- 
ing gap, the reduction of the stiffness at finite tempera- 
ture is sharper than in the uniform superconductor. It 
therefore does not immediately follow that the interplay 
of these two effects can lead to an enhancement of the 
critical temperature. Previously we found that such an 
enhancement is possible in a bilayer of underdoped and 
overdoped cuprates, under appropriate conditions"'. In 
the present scenario, however, we find that of the inho- 
mogeneously doped layer is lower than the one of a homo- 
geneous underdoped superconductor of doping p\. The 
reason is that already at T = the enhancement of the 
stiffness due to enlarged carrier density is counteracted to 
a large extent by a significant paramagnetic suppression 
of the stiffness which is inevitable in inhomogeneous su- 
perconductors. Consequently, the zero temperature stiff- 
ness is enhanced compared with the uniform case, but 
not enough to allow for an enhancement of Tg. 

Nonetheless, we find that the critical temperature of 
the system increases with the reduction of the relative 
width of the metallic stripes. This allows for a large 
proximity gap in the metallic regions, manifested in a 
relatively small reduction of the stiffness at finite temper- 
ature. In order to maximize the proximity effect, but at 
the same time allow for a significant contribution of car- 
riers from the metallic region, an optimal configuration 
should have small but relatively dense metallic regions. 
In the following we consider the effect of small metallic 
regions. 
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B. Microscopic Inhomogeneity 

111 this model, described in Sec. IV we assume micro- 
scopic overdoped regions (doping P2) which are placed in 
a low doping superconducting background (doping pi), 
see Fig. 1(b). The microscopic overdoped regions are 
modeled as single site impurities with zero or very weak 
local Hubbard repulsion {U ~ 0), which induces modified 
hopping and pairing amplitudes along their neighboring 
bonds, as depicted in Fig. 4. The hopping amplitude 
along these bonds is the bare t rather than the renor- 
malized value of SBMFT, and the local pairing strength 
there is suppressed to zero. 

In the presence of the bond disorder we compute the 
temperature dependent superfluid stiffness using a per- 
turbative expansion to second order in the impurity 
strength for disorder averaging (second order Born ap- 
proximation). Then we determine the transition temper- 
ature using the Kosterlitz-Thouless criterion as before. 

Since the variations in doping generates unconven- 
tional bond disorder, the calculation bears several impor- 
tant differences from the standard impurity averaging. 
The most important difference is that the bond disorder 
introduces local modulations in the current operator thus 
renormalizing the coupling to the external vector poten- 
tial. As a result, the superfluid response obtains vertex 
corrections which have no counterpart in standard (on- 
site) impurity averaging but play a crucial role in our 
case. One important effect of these corrections is to allow 
for an enhancement of the zero temperature diamagnetic 
stiffness of the disordered system compared with the pure 
one. A second effect of the vertex corrections is to intro- 
duce a paramagnetic reduction of the stiffness at zero 
temperature, similarly to the mesoscopic inhomogeneous 
scenario. In addition, the disorder introduces self-energy 
corrections which amount to an anti-proximity effect that 
acts to reduce the average pairing gap and contributes to 
the suppression of the stiffness at finite temperature. 

The net effect that we find is an enhancement of the 
superfluid stiffness and a concomitant increase in the crit- 
ical temperature for a given bulk doping level pi. In- 
terestingly, we even find an overall enhancement of the 
maximal Tc, that is at optimal doping, compared to the 
maximal Tc of the homogeneous system. 

III. MESOSCOPIC SCALE INHOMOGENEITY 

A. The Model 

In this section we consider a stripe model. The inho- 
mogeneity is of mesoscopic scale in the sense that the 
width of the stripes is of the order or somewhat larger 
than the coherence length associated with the supercon- 
ducting regions. The superfluid response of such a striped 
system is of course anisotropic. However we envision that 
it becomes isotropic on macroscopic scales due to mixing 
of striped domains with random orientations as sketched 



in Fig 2(a). The doping level of the stripes alternates 
between pi in underdoped superconducting (SC) stripes 
of width I, and p2 in metallic (M) stripes of width d. 

As the Hamiltonian of a single domain we take the t— J 
model 

H = PG^%4cjaPG+5I>^uS,S, (3) 

where Pq is the Gutzwillcr projection that eliminates 
double occupancy of sites in the superconducting stripes, 
but does not affect the metallic stripes. The magnetic ex- 
change coupling is Jij = J in the superconducting stripes 
and it vanishes in the metallic stripes. 

We treat the space dependent projection and ex- 
change interaction within slave boson mean field the- 
ory (SBMFT)"'''-''. The resulting Hamiltonian is of the 
form (2), with space dependent Xij ^^'^ ^ij- The 

electro-chemical potential Hi is determined such that the 
doping levels of the superconducting and the metallic re- 
gions are pi and p2 respectively. Due to the spatial varia- 
tions in doping the renormalization of the hopping varies 
in space too and equals g\j = 2pi/{pi + 1) in the super- 
conducting stripes and gl^ = 1 in the metallic stripes, 
while the tunneling at the interface between the two re- 
gions is renormalized by gl^ = \j2p\l{p\ -I- 1)). 

Given all the parameters of the mean field model, the 
fields Xij s-iid Aij can now be determined by the self 
consistency conditions: 

3 \ , ^ . 

CJ 

A. = ^(^It4i-i4t) (4) 

An example of the resulting profile of the pairing ampli- 
tudes is plotted in Fig. 2(b), where A^^ and Aj, denote 
the pairing amplitudes on bonds along the x and the y 
directions respectively. Because the pairing amplitude in 
the metallic regions is non zero, these regions contribute 
to the superfluid stiffness at low temperatures. 

B. Calculation of the Superfluid Stiffness 

In a striped system the superfluid response depends 
on the direction of the applied phase twist. However, 
we assume that the system consists of many striped do- 
mains with random orientations along the principal axes. 
Under this assumption the superfluid response is homo- 
geneous on large scales. It was shown in Ref. 32 that 
the superfluid stiffness of the mixed domains is given by 
the geometric mean of the x and y components of the 
stiffness of a single domain ps = \J KxxKyy. Here Kaa 
(a — X, y) are the diagonal components of the response 
tensor, 

A-..^^ (5) 
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FIG. 2: (a) A model of the inhomogeneous layer as an array of striped domains, which on average is macroscopically fourfold 
rotationally invariant, (b) The self consistent gap profile, solved for J — t/3,t' — 0,x = 0.1, d — 3a. SC and M denote the 
underdoped and overdoped regions respectively. 



where A6'f, is the static phase difference applied across 
the system in the h direction and la = J J ads a is the 
total current measured in the a direction. 

In an inhomogeneous system, we express the stiffness 
tensor using the microscopic response kernel K,ab{fi f') de- 
fined through 



J air) 



dr' Kab{r,r')dbe{r') 



(6) 



t<t' 



The response kernel can then be calculated using the 
standard Kubo formalism. In the y direction, parallel 
to the stripes, the stiffness Kyy — J dxJy/ AOy is sim- 
ply an algebraic sum of the response kernels along the x 
direction 



Kyy = 



dX I dX tvyyi^X ^ X , Cj^y O); 



(7) 



where we used the uniformity along the y direction to 
express it in terms of the gj, = Fourier component of 
the response kernel. 

To derive an analogous relation between Kxx and 
Kxx{x, x'), we follow the arguments presented in Rcf. '-VI. 
It is convenient to use the lattice formulation and ex- 
press all convolution integrals as matrix products. The 
response kernel is then defined by (summation over re- 
peating indices implied), 



Jaii) = Kab{ij)db9{j) 



(8) 



where the y dependence is suppressed and we denote by 
i,j the position in the x direction only. 

A static current is divcrgcnceless V • J = and there- 
fore derived from a potential, Ja{i) = eabdb<Pii)- Plugging 
this back into (8) we obtain 



iabdb(t>{i) = Kab{ij)db0{j). 



(9) 



We can now derive a second relation between and 9 
if we multiply by eK~ 
iacKZledb we arrive at 



from the left. Defining 



I^abij^)9b4'{i) = eabdbO{j). 



(10) 



Equations (9) and (10) establish a duality relation 9 ^ 
<j> and K o K^. We make use of this duality in the 
calculation of Kxx- The response Kxx is obtained by 
applying a phase difference A^^; along the x direction 
and measure the resulting current Jx, equivalent to a 
difference in (j> in the transverse direction, Acjiy = J dyJx- 
The response is then Kxx = A(j>y/A9x 



Using Eqn. (10) we deduce that dx9 



^yy^y 



dy(j}, 



which allows us to apply relation (7) with 9 (j) ^-nd 
Doing so, we obtain the response A9x/A<f)y = 



K K 



D 



J dxdx' K. 



D 

yy 



Kx 



This gives the result 



J dx J dx'Kxx {x,x') 



(11) 



When the stripes are macroscopic wc can take the re- 
sponse functions to be translationally invariant within a 
stripe. Then (11) reduces to the well known fact that the 
stiffness of macroscopic objects in series adds like resis- 
tors in parallel. 

The Superfiuid stiffness is now expressed in terms of 
the local response Kernel which can be computed using 
the standard Kubo formalism. The diamagnetic contri- 
bution of the response to a transverse vector potential 
is 

Kf^^{x,x',qy^O) = eHxx'{-Ka{x,qy^O)). (12) 

where Ka{j) = -t{j) E<t(c]+x„,^Cj- ^ + 
local kinetic energy operator. The paramagnetic contri- 
bution is 



(13) 



i\,b {x,x ,qy,iu}n) = dre*""^ X 

Jo 

X {T,jP{x,qy,T)jP{x',-qy,0)), 

in the limit lu = 0, gj,— >0. Here is the paramagnetic 
current operator: 



J 



(14) 
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In order to calculate (12) and (13) we diagonalizc the 
Hamiltonian using the Bogoliubov transformation 

a 
a 

We then solve the self consistent equations (4) and ex- 
press Ka{j) and in the new basis. As an example, the 
X component of the static paramagnetic response kernel 
is 

kyaP 

Ea.ky ~ Ep^ky—qy 

where V denotes the principal part, Ea^ky are the eigen- 
energies, /(e) = [1 + e^*^]"-^ is the Fermi-Dirac distribu- 
tion, and 

Jif{j,k,q) = tt{j)[ul{j + l,k)up{j,k-q) 
+ <(j + 1> k)vp{j, k~q)]+ c.c. 

To solve the periodic problem we introduce an additional 
superlattice momentum, whose index is suppressed here 
for simplicity. 

Finally, the superfluid stiffness Ps(T) = ^ K^^Kyy 
is computed from the response kernels Haaix^x') = 
nf^ix.x') + kPI'^{x,x') using Eqns. (7) and (11). 



C. Results and Discussion 

We first discuss the superfluid stiffness at zero temper- 
ature, and then turn to an analysis of the temperature 
dependence of the stiffness, in order to estimate the crit- 
ical temperature of the striped system. 

The superfluid stiffness at zero temperature Ps(T = 0) 
is plotted in Fig. 3(b) as function of the relative width of 
the metallic stripes, C = d/{d + I) at flxed width I + d = 
20a {a is the lattice spacing). The doping levels of the 
superconducting and the metallic stripes are pi = 0.1 
and p2 = 0.35 respectively. Note the enhancement of the 
stiffness compared with the uniform superconductor pgc 
which is maximized for ^ 0.5. 

The enhancement of the zero temperature stiffness and 
the optimal volume fraction arc determined by the in- 
terplay of two effects. First, the metallic regions are 
gapped by the proximity effect, and contribute their large 
number of carriers to the diamagnetic superfluid den- 
sity K^l"" which increases with However, this increase 
is partially countered by a zero temperature paramag- 
netic term K^y'^'^ special to inhomogeneous supercon- 
ductors. A similar effect was noted by us in a bilayer 
het erostructure ^ ' ' . 



To see if the moderate net increase of the zero tem- 
perature stiffness will facilitate enhancement of the tran- 
sition temperature we compute the full temperature de- 
pendence of the stiffness. As an example Fig. 3(a) shows 
the result for a specific ratio d/{d + I) = 3/20 with 
Pi ~ 0.1 and p2 = 0.35. The macroscopic stiffness ps 
is seen to decrease approximately linearly with tempera- 
ture, as in a uniform d-wave superconductor, but with a 
larger slope dps/dT. As a result, the transition tempera- 
ture, determined using the Kosterlitz-Thouless criterion 
Ps{Tc) = 2Tc/tt, is found to be lower in the inhomoge- 
neous layer despite the increased stiffness at zero tem- 
perature. This remains the case in all possible stripe 
geometries, as shown in Fig. 3(c). 

The slope dps/dT is affected by two main factors: the 
first is the density of states (DOS) of low energy quasi- 
particles that carry the paramagnetic current and the 
second is the effective charge of these quasiparticle (or 
the effective current renormalization). The DOS of the 
system is plotted in Fig. 3(d). Below a threshold energy 
oi E ^ 0.05t the DOS is the same as in the uniform super- 
conductor of pi =0.1. This agrees with experimental re- 
sults in inhomogeneous cuprate superconductors '. Note 
that the limit of very narrow metallic stripes (d = 3a) 
preserves the low energy DOS up to a relatively high en- 
ergy, compared with the critical temperature. As the 
metallic stripes get wider, there are more low energy 
states that contribute to the reduction of the stiffness 
with the temperature. 

Despite the fact that the low energy DOS is the same as 
in the uniform superconductor, the slope \dps/dT\ is still 
larger than in the uniform case. This is a consequence 
of the difference in the effective charge of quasiparticles 
in the two systems: In the underdoped superconducting 
regions the quasiparticle charge is renormalized down by 
a factor of gt = 2pi/{\ +pi), whereas in the metallic 
regions there is no such renormalization and the current 
is carried by electrons. As a result, at finite temperature 
the stiffness reduction in the inhomogeneous system is 
steeper than in the uniform underdoped superconductor. 

Here we should note again that, in general, the renor- 
malization of the current carried by a quasi-particle that 
enters the low temperature dependence of the stiffness 
is a Fermi-liquid parameter that may be renormalized 
compared to the SBMFT value of = 2pi/(l + pi). In- 
deed measurements of the temperature dependent stiff- 
ness give a renormalization factor is independent of dop- 
ing over a wide range of doping in contradiction to the 
mean field prediction. However for an inhomogeneous 
system there is no unambiguous way to replace the mean 
field value of the current renormalization by a single phe- 
nomenological parameter. Moreover the existence of a 
length scale d of the superconducting regions may intro- 
duce a cutoff that prevents this parameter from flowing 
far from its mean field value at low energies. 

The stripes model shows that doping inhomogeneity 
on a mesoscopic scale can lead to an increased superfluid 
stiffness at zero temperature. This is a consequence of 
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FIG. 3: (a) Temperature dependent superfluid stiflFness of a striped versus a homogeneous superconducting layer. In this 
example the width of metallic stripes (doping p2 = 0.35) in this example is d = 3 and that of the superconducting stripes 
(J = t/3 and doping pi = 0.1) is / = 17 (in lattice constants). The critical temperature estimated from the Kosterlitz-Thouless 
criterion, is seen to be higher in the homogeneous system (circle marks), (b) Zero temperature stiffness vs the relative size 
of the metallic segments, d/L for a fixed L = 20. The diamagnetic (dashed) and paramagnetic (dot-dashed) contributions to 
the stiffness Kyy axe plotted for comparison. All results are normalized by the zero temperature stiffness psc of the uniform 
superconductor, (c) Critical temperature, normalized by the uniform value, vs d/L for the same parameters as in (b) and 
different unit cells L (in units of the lattice spacing a), (d) The density of states for a system with the same parameters as 
(a). At low energies the density of states (DOS) of the striped system (black) is identical to the DOS of a uniform underdoped 
superconductor of the same size (gray). The dashed line denotes the calculated slope dv/dE aX E = Q for the uniform 
superconductor. 



a proximity gap that opens in the metallic stripes which 
then contribute their high carrier density to the stiffness. 
On the other hand the metallic stripes also give rise to 
low energy states that hasten the reduction of stiffness 
with temperature. In addition there is a paramagnetic 
reduction of the stiffness even at zero temperature due to 
the impurities. For these reasons the transition temper- 
ature of the striped system is found to be always lower 
than that of the homogeneous system. The highest Tc is 
obtained for the narrowest metallic stripes because then 
the proximity coupling to the bulk is high and the An- 
dreev bound state are only slightly below the gap. It is 
therefore tempting to consider the case of even smaller 
metallic regions by reducing the length of the stripes in 
addition to their width to a microscopic scale. In the 
following section we examine a model that takes a step 
in this direction. 

IV. MICROSCOPIC SCALE INHOMOGENEITY 

A. The Model 

In this section we consider a scenario in which the 
metallic regions embedded in the underdoped supercon- 
ductor are nearly point like. We model these metallic 
impurities as cross vertices of the square lattice (see Fig. 
4) on which the average doping p2 is higher than the 



bulk average pi . The effective hopping matrix elements 
and the pairing amplitudes on these links naturally take 
different values than the bulk. Specifically, in the mean 
field model of Eqn. (2) the parameters gjj, Aij and Xij 
take a different value on the impurity bonds. 

We analyze two scenarios: in the first, the impurities 
are metallic, with doping p2 > 0.3, such that on the 
impurity bonds gij = 1 and Aij = Xij = 0. In the other 
scenario the excess doping on the impurity sites P2 —pi is 
small, leading to gij ~ g{p2) ^^nd Aij ~ A{p2) with the 
doping dependence of SBMFT. In this case we assume 
that X, which has a very weak doping dependence, is 
uniform throughout the system. 

The Hamiltonian H = Hq + Himp consists of a uni- 
form part and an impurity contribution. Written in mo- 
mentum space, the uniform Hamiltonian is the Fourier 
transform of (2), 

Ho = Yl *k(?ka3 + Au(Ti)*k. (17) 

k 

Here ^'j^. = {'4ct' ^-kj.}, are Pauli matrices, = 
— 2icff (cos kx + cos ky) — p, and Ak = Ao(cos k^^, — cos ky) 
with icff = g{pi)t + X- The impurity Hamiltonian is 

i^imp = *kPk'-k [Uu,'<7z + Vkk'CTl] *k' (18) 
kk' 
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B. Calculation of the Superfluid Stiffness 

The stiffness is the linear response of the system to 
an externahy applied vector potential A(r). In order 
to calculate it in the disordered system, it is conve- 
nient to resort to the real space Hamiltonian (2) and 
include a vector potential through a Peierls substitution, 
gijt^gijtcxp[ieAij] = gijtcxp[ieAx{ri)] in the case of 
a vector potential in the x direction. For the linear re- 
sponse calculation we expand the Hamiltonian to second 
order in At'''', 



FIG. 4: Doping inhomogeneity on a microscopic scale. An 
illustration of the model. Nearest neighbouring bonds to the 
impurity sites (solid lines) are characterized by enhanced hop- 
ping amplitude gt + St and reduced pairing Ao -I- 5a, with 
respect to the superconducting background of gt and Aq re- 
spectively. 



where 



Ukk' = — 2(5t(cos kx + cos ky + cos k'^ + cos k'^) 
Vkk' = Sa (cos kx — cos ky + cos k'^ — cos ky) 
Pk'-k = ^e-^O^'-'')--" (19) 



Here Sa = A(p2) — A(pi), and {va} are the impurity 
sites. The excess hopping at the impurity sites is 6t = 
tcs{p2) — icff(pi)- In the case of metallic impurities we 
set ioff(p2) = t. 

The above terms result from the shift in doping level 
from pi to p2 near the impurity. We should in principle 
include also the direct impurity potential, which caused 
the change in hole concentration. Such a potential that 
acts locally on the impurity as U{r) = Uo5{r — r^), 
can be regarded as a fc independent contribution to 
C^kk'- The magnitude of this term can be estimated 
from the observed change in hole concentration through 
Uq — {p2 — Pi)/^, where k is the local compressibility. 
We omit this term from the calculations described below. 
Then, at the end of Sec. IV C we quantify the contribu- 
tion of the direct impurity potential and explain why it 
can be neglected. 

Our goal is to compute the temperature dependent su- 
perfluid stiffness and estimate the transition temperature 
of the inhomogeneous layer compared to a uniform layer. 
To this end we use the Born approximation to perform 
the disorder average. This is strictly valid in the limit of 
dilute uncorrelated impurities and weak disorder. We ex- 
pand to first order in the impurity concentration Ui and 
second order in the strength of a single impurity St/tos 
and (5a/ Aq. In practice we will allow Sa/^o to be close 
to —1 which is the case when the overdoped inclusions 
are already in or close to the metallic regime. 



HiA,) = HiO) - J2 



eMr)A,{r) + -lUr)Alir) 



(20) 



with 



jx{r) = i^ta:(r)(4+^ ,,Cr,^ - 4_^Cr+x,^) 
r,(T 

Kx{r) = -^t:,(r)(4+^„Cr,<x +4,^Cr+x,<T)- 



Here tx{r) = g{pi)t + J2a ^("^ ~ ^a) is the coupling to 
the external vector potential, in the presence of the mod- 
ified bonds around sites r^. The excess local current on 
the impurity sites is = t[g(p2) — g{pi)]. In the case of 
highly overdoped impurities (p2 > 0.3) we take 5(^2) = 1- 
Note that this impurity contribution is different from St 
that appears in the impurity Hamiltonian (18). The rea- 
son is that the external vector potential couples only to 
the hopping gijt, and not to the Fock term proportional 
to Xj which originated from the magnetic exchange inter- 
action. 

The superfluid stiffness is now given by 



Ps 



lim Uxx{ci,uj = 0) 

q— >0 



(21) 



where X denotes the average over disorder realizations 
and. 



n:cr(q,w„) 



/3 



dre*"""(jx(q,r)j,(-q',0)) 



Note that after disorder averaging the RHS is propor- 
tional to <5qq'. The different contributions to ps are pre- 
sented as diagrams in Fig. 5, where we denote diamag- 
netic terms by Da and paramagnetic terms by Hq . 

One type of correction to the stiffness stems from stan- 
dard rcnormalization of the electron self-energy by the 
impurities. Such corrections are given by diagrams Dq 
and Hq in Fig. 5. Similar terms would arise in the com- 
mon case of point (on-site) impurities. We note that the 
vertex correction Hi vanishes due to inversion symmetry 
of the impurity potential. 

A second type of correction to the stiffness is special to 
the bond disorder we consider here. The disorder in the 
hopping amplitude introduces modulations in the local 
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FIG. 5: Diagrammatic calculation of the superfluid stiffness in 
the disordered system. Top: diamagnetic contributions. Mid- 
dle: Paramagnetic contributions. Bottom: Green's function 
renormalization and the definition of the self energy within 
Born approximation. A dashed line corresponds to a scat- 
tering event and x denotes a single impurity. Note that the 
scattering is a matrix in Nambu space. 



current operator and kinetic energy, proportional to St- 
This causes a direct renormalization of the coupling to 
the external vector potential, as represented in diagrams 
£'i,L>2,n2,n3 and in Fig. 5. 

Self energy corrections. The disorder in the hop- 
ping and pairing strength introduces renormalizations to 
the spectrum parameters or to the electronic Green's 
function, which in turn, affect the superfluid stiffness. 
Such corrections are represented by diagrams Dq and 
flo in Fig. 5. In order to calculate these diagrams we 
first compute the rcnormalized Green's function using 
the Born approximation. 

The Dyson equation for the disorder averaged Green's 
function is depicted in Fig. 5 and given by 



where the bare Green's function is 



(G'k,tj„) 



(22) 



(23) 



and Sk,Lj„ = X]a=o '^aCa is the self-energy after disor- 
der averaging. To calculate the self energy explicitly 
in the limit of small w„, we use the fact that the main 
contributions arise from the vicinity of the nodal points 
Ak = ^k = 0. Wc expand around these points and solve 
self consistently for the decay rate Eg in the limit of 
w— >-0, k— )-knodc similarly to Ref. >)4. The high energy 



cutoff for this approximation is defined as po- This cal- 
culation gives (see Appendix A) , 



So 
Si(k) 

S3(k) 



2n,^(l-77) Ak 
2ni — (1 - 7^)^k + Sf, 



(24) 



where 



V 



S^Sa/^o + (1 - Si)6t/t,s, Si 



■\/2Ao/ioff/7r = '2 \/vA/vf / IT , and (5^ is a k independent 
constant that renormalizes the chemical potential. 

The low energy limit of Eq is exponentially small close 
to zero doping, and is further suppressed by the large 
number tlg/{ni5f). We solve for the other components 
of the self energy under the self consistent assumption 
that any k dependence of Sq is negligible and indeed get 
that the entire effect of the decay rate iSq is negligible. 
For more details about the calculation the reader should 
turn to Appendix A. 

In the absence of decay, no zero energy states are in- 
troduced by the disorder. The effect of Si and S3 is to 
rcnormalizc the gap and the hopping, leading to a cor- 
rected spectrum i?k = [^k + ^kl^^^- energy 
limit this is equivalent to a renormalization of the effec- 
tive values oi Vf and va which we find to be, 



Vf = Vf 



VA = VA 



1 + 2n,-^ (I - 77) 

toff 

l + 2n,^(l-,7) 
^0 



(25) 



The renormalization of va is the anti-proximity effect due 
to the metallic inclusions, which gives rise to a modified 
coefficient of the linear DOS compared with the pure sys- 
tem. These modifications primarily affect the low tem- 
perature physics in the disordered system. 

With the Green's function at hand we can calculate the 
leading contributions to the superfluid stiffness. Details 
of the calculations appear in Appendices B and C. The 
contributions to the superfluid stiffness, to second order 
in the disorder strength, can be separated into zero tem- 
perature and flnite temperature contributions. 

Zero temperature.- The contribution to the zero tem- 
perature stiffness due to self energy corrections is the 
diamagnetic response expressed in diagram Dq. This is a 
non-universal contribution which turns out to differ only 
very slightly from the bare diamagnetic stiffness of the 
pure system (see Appendix B for details). 



Dq = 2gt cos kx 



1 - 



6 



gtVo [ ^ 



(26) 



where VqIX) is an order unity slowly decreasing function 
of its argument in the relevant range of parameters. 



10 



Note that in practice, this function may include a weak 
dependence on the chemical potential which we neglect, 
assuming low doping. To conclude, the renormalization 
of the spectrum parameters due to the self energy correc- 
tions have a negligible effect on the diamagnetic stiffness. 

Finite temperature. - The finite temperature contribu- 
tion to the stiffness due to self energy corrections arises 
from diagram Ho. The effect of disorder here is to modify 
the low energy density of states through a renormaliza- 
tion of the effective values oi Vf and v^- This affects the 
superfluid stiffness through the paramagnetic contribu- 
tion fig leading to faster reduction of the stiffness with 
temperature. More precisely 



2 log 2 8(Zi)2 



VfVA 

Ao 



T 



(27) 



^cff 



boT. 



Here ~ —dpg/dT is the slope in the clean system and 
Z is the renormalization of the quasiparticle current by 
interactions. The low T behavior is dominated by low 
energy quasiparticles, which may be altered by Fermi- 
liquid renormalization not included in the mean field the- 
ory. Therefore Z should be taken as a phenomenological 
Fermi- liquid parameter^''^'*'^ and not as the value g{p) dic- 
tated by the microscopic mean field theory. 

In our case the disorder acts to induce faster decrease 
of the superfluid stiffness with temperature. This is 
because when the inclusions are highly overdoped with 
nearly zero gap then Ja/Aq ^ — 1, while < St/tcff ^ 1. 

Current operator renormalizations. The second 
type of corrections to the stiffness have no counterpart 
in systems with standard on-site disorder. The disorder 
in the hopping amplitude introduces renormalizations 
of the kinetic energy and the current operator 
jx, proportional to S^. This leads to corrections of 
0{S^),0[{S^)'^] to the stiffness, that are represented 
as vertex corrections in diagrams _Di, _D2, 112, Ha, 1X4 of 
Fig. 5. We again distinguish between zero tempera- 
ture and finite temperature contributions to the stiffness. 

Zero temperature.- The most intuitive effect of the 
vertex correction is the increase of the diamagnetic stiff- 
ness at the impurity sites due to the extra charge carriers 
they contribute. This effect is reflected in the diagram 
Di with each impurity bringing an additional 26^ to the 
average kinetic energy 



Di = 4niSf ^ cos fcx 1 1 — 
k V 



Vf 



gt 



(28) 



This expression reveals a small parameter, UiS^/gt, that 
did not appear in the Hamiltonian. The perturbative 
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FIG. 6: Superfluid stiffness and critical temperature in the 
inhomogeneously doped layer with point impurities. In pan- 
els (a) and (b) the impurities liave an average doping charge 
P2 = 0.35 while the bulk doping pi varies, with impurity con- 
centration Ui = 0.1. (a) Zero temperature stiffness (Solid) 
compared to that of the clean case (dashed). The Dashed- 
dotted line marks the diamagnetic contribution, (b) Tc with 
impurities compared to without. We used a doping in- 
dependent quasiparticle current renormalization of ^ = 0.5. 
Results are plotted only within the validity range of the dia- 
grammatic expansion, (c) Relative change in Tc for the case 
of small excess doping on the impurity, with impurity concen- 
tration Hi = 0.2. Contours map the relative change T^/T^ as 
a function of the base doping pi and the excess doping on the 
impurities p2 — pi . 



correction inevitably becomes large upon underdoping 
towards the Mott insulator where riiS^/gt — > 00. This 
signals the breakdown of the Born approximation at dop- 
ing levels smaller than pl ~ g{p2)ni/{l + 2ni). 

The second significant contribution to the zero tem- 
perature stiffness stems from the paramagnetic diagram 
114, which is seen to be 



Hi = -2n,V^^o 



off 



VA 
Vf 



(29) 



Here V^iX) is an order unity decreasing function of its 
argument in the relevant parameter range. This term is 
closely analogous to the zero temperature paramagnetic 
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reduction in the stripe model of sec. III. Here as in the 
stripe model, The effect acts to moderate the enhance- 
ment of the stiffness at zero temperature. 

Another correction to the zero temperature stiffness is 
given by the diagram D2. This diagram, which repre- 
sents a combined renormalization of the vertex and the 
spectrum, is calculated to be 

[Ao \vf J toff \vf J 

Here 2?i(X) is an increasing function and 'D2{X) is slowly 
decreasing, and their weak dependence on the chemical 
potential is again neglected. This diagram turns out to 
give a negligible numerical contribution to the overall 
stiffness. 

Finite temperature. - The finite temperature contribu- 
tions to the stiffness that arise from current renormal- 
ization are shown in diagrams H2 and H3. An explicit 
calculation gives 

H, . -4n,^^ T 

gt TT W/WA 

2 log 2 {2V2Ztf ^ 

Hs = Arii 77 — T. (30) 

gt IT VfVA 

Within SBMFT the current renormalization Z ~ g{p) 
depends strongly on the doping. However, it is known 
that this strong doping dependence leads to a disagree- 
ment with the experimentally measured slope dps/dT, 
which is seen to be almost independent of doping."" 

Here we adopt a phenomenological approach, with an 
effective paramagnetic current renormalization Z which 
is independent of doping.-'^''*" This holds at finite low 
temperature, when the physics is dominated by the ef- 
fective theory of low energy Dirac quasiparticles. In this 
case, the entire contribution H2 + II3 is negligible be- 
cause it stems precisely from the difference in the local 
current operator between the pi superconductor and the 
P2 impurities. Therefore, when summing up the finite T 
contributions to the stiffness we neglect these two dia- 
grams. 

C. Results and Discussion 

We can summarize the results of this section by putting 
together the various contributions to the superfluid stiff- 
ness. This gives the temperature dependent stiffness 

p,(T) = po - boT + 2n,(Jp,(0) - Sb T) (31) 

Here the first two terms constitute the usual expression 
for the temperature dependent superfluid stiffness of a 
uniform d-wave superconductor-^'', as reviewed in sec. II. 
The second term is 

^^^(O)-^Po-^^o(^-j (32) 



The leading order correction to 5ps{Q) in the impurity 
strength is due to the added charge carriers donated 
by the impurities. The negative second order term is 
a paramagnetic correction to the zero temperature stiff- 
ness analogous to the paramagnetic correction that we 
derived previously for a bilayer heterostructure. In the 
latter case this correction was proportional to (Ji~- J2Y' , 
the square of the difference of the quasi-particle currents 
on the two layers. Here similarly this contribution scales 
as ((5()^ a [g{P2) — g{Pi)\^i which is the square of the dif- 
ference between the local current renormalization in the 
bulk and near the impurity. 

The last term in Eq. 31 is the change of the linear 
reduction of the stiffness with temperature due to the 
presence of impurities. It is given by 

where Si = 2-\/?jA/w//7r and &o is the parameter for 
the uniform superconductor given in sec. II 60 = 
2\og2{2y/2Zt)'^ / {-KV fV/^Y'' . Note that the expression in 
the first bracket is positive because < on the impu- 
rities. Hence the superfluid stiffness is reduced faster as 
a function of temperature than in the uniform supercon- 
ductor. We estimate the parameters of the uniform sys- 
tem using SBMFT, so that Vf ^ 2^/2^off and = V2A. 
Taking J = t/3, the effective hopping and gap parame- 
ters are given by toff(p) = gip)t + X ^-nd A(p) = x[l — 4p], 
where x is the value of the mean fields (both pairing and 
Fock field) at zero doping. 

Figure 6(a) displays the calculated zero temperature 
stiffness as function of the doping pi with the impurities 
fixed to a high doping level p2 = 0.35, which corresponds 
to zero pairing amplitude, and a hopping amplitude of 
t. We plot the total stiffness PsiT = 0) as well as the 
diamagnetic contribution pdia — Dq + Di + D2. Note that 
the diamagnetic contribution in the disordered system 
Pdia is significantly increased with respect to the pure 
case, pq. However, the total zero temperature stiffness 
Ps{T = 0) is only moderately increased compared to the 
uniform case (where po is the total stiffness at T = 0). 
The reason for this is the zero temperature paramagnetic 
contribution of the impurities H4 . 

In panel (b) of the same figure we plot the critical 
temperature as a function of the bulk doping pi, esti- 
mated from Eq. (33) using the Kosterlitz-Thouless cri- 
terion Ps{Tc) = 2Tc/n. Again this is for a fixed value 
of P2 ~ 0.35 and n,j =0.1 and the result is compared 
against of the pure system. The critical temperature 
of the disordered system is significantly enhanced, above 
the maximal Tc of the clean superconductor. The max- 
imum of Tc is shifted to the underdoped regime. These 
results are reminiscent of experiments by Yuli et al' ' that 
show a Tc enhancement in an underdoped-overdoped bi- 
layer. We can relate our results to the experiment if we 
assume that the interface between the two layers is in 
fact an inhomogeneous mixture of underdoped and over- 
doped regions. Our results suggest that an optimal con- 
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figuration for Tc enhancement can be achieved by placing 
point-like metallic inclusions inside a slightly underdoped 
superconductor . 

Figure 6(c) shows the relative change in the critical 
temperature with respect to as function of pi and 
Sp ~ P2 ^ Pi, for Jii ~ 0.2. The critical temperature is 
enhanced relative to the clean system by up to 15%, in 
a broad range of pi and Sp. Here the excess doping on 
the impurities Sp is small, and there is no enhancement 
of Tc above the maximal Tc of the clean superconductor. 
The main reason for this is the zero temperature param- 
agnetic reduction of the stiffness due to the impurities. 
Without this effect we could have obtained an absolute 
enhancement of Tc in the disordered system, also in the 
small Sp limit. We have checked and found that whether 
we use the microscopic or phenomenological parameter 
Z to renormalize the quasiparticle current makes very 
little difference to the final result of Tc- 

It is instructive to look at the behavior of the stiffness 
and Tc, for small values of Sp = p2—pi, for which we can 
neglect second order contributions in St/t^s and |(Sa/^oIj 
such as the paramagnetic effect. Here we use the SBMFT 
doping dependence for both the bulk and the impurity, 
such that St ~ Si = t[g{p2) — g{pi)\- In this regime there 
is a simple expression for the superfluid stiffness. 



Ps{T) ~ Pc\{T) + 2ni-^po - 4ri 

gt 



Sa 
Ao 



^]boT 

toff , 



where Pc\{T) ^ pa — boT is the stiffness of the clean su- 
perconductor. The zero temperature stiffness is always 
enhanced, whereas the slope \dps/dT\ is increased. The 
latter is easily seen by expressing |(5a/Ao| — ^t/ieff as 
function of pi and Sp. Using the Kosterlitz-Thouless cri- 
terion as above we can estimate the change in transition 
temperature Tc with respect to the critical temperature 
Tc of the clean superconductor. 



TP 



1 + 2n,; 



gt 



2bo f St 


Sa 


)] 


6o + 2/7r \t^s 


Ao 





(33) 



This can be expressed in terms of the doping level pi 
of the clean superconductor and the difference in doping 
Sp = P2—pi between the background and the impurities. 
We obtain an expression of the form 

l£ = l + 2r.,^(p*-pi)J-(pi). 
T^c Pi 

where J^{pi) is positive for pi < 0.25. This implies 
that for pi < p*, the critical temperature of the disor- 
dered system is enhanced compared with the clean su- 
perconductor with doping pi . Under the assumptions of 
SBMFT, with Z = 0.5 and x = OAt, we get p* - 0.125 
and 



(l-pi)(0.08+pi)(0.6+pi 



(0.5-pi)(0.25-pi)(0.17+pi)(0.29+pi)' 



(34) 



We note that J'(O.l) ~ 18. 

It seems from Eq. (34) that Tc can be further 
enhanced by increasing the impurity concentration. 
However, by doing this we would quickly violate the 
Born approximation. In particular, the superfluid stiff- 
ness in this non perturbative regime should be calculated 
as the resistance of an effective resistor network with pj^ 
of the various puddles playing the role of the resistance. 

We now remark on the nature of our perturbative ap- 
proach and the small parameters involved in it. The scat- 
tering from individual impurities is taken into account 
within the Born approximation to second order in the 
impurity strength as measured by the parameters St/teff 
and |(5a/Ao|. We found that the second order correction 
to both ps and T^ was always negligible compared to the 
first order contribution. This was the case even when we 
took the parameter |(5A/A| ~ 1. An additional small pa- 
rameter riiS^/gt appeared through the effect of the impu- 
rities on the coupling to the electromagnetic field rather 
than the scattering on the impurity potential. 

We would like to contrast our approach with the 
commonly used self consistent T-matrix approximation 
(SCTMA) '^', which treats the single impurities exactly. 
This turns out to be important to describe the effect of 
in-plane ion substitutions such as Zn impurities that act 
as unitary scatterers and give rise to strong bound states. 
However in our case the SCTMA is not analytically solv- 
able because of the strong momentum dependence of the 
impurity potential and the fact that it acts as a matrix 
in Nambu space {St is the diagonal component and 
off-diagonal). Fortunately, the disorder potential that 
interests us is generated by dopants, which reside out- 
side the CuO plane. Indeed we can show that the fact 
that such impurities induces only a small change in the 
local doping level {p2—Pi ^ 1) implies that the impurity 
scattering is far from the unitary limit and does not give 
rise to a bound state. To see this consider the strength of 
the local impurity potential Uq ~ {p2 —pVj/n. Since the 
compressibility k is approximately the density of states 
at the Fermi level fo , the dimcnsionless impurity strength 
is just UqVq ~ P2 — pl ^ 1- This is far from satisfying 
the condition for formation of a bound state. In this 
limit the direct impurity scattering can be taken within 
the Born approximation and lead to negligible contribu- 
tions to the low energy DOS. Hence it leads to a 
concomitantly small correction to the stiffness. 



V. CONCLUSIONS 

We investigated the effects of doping inhomogeneity 
on the superconducting properties of the cuprates using 
the slave boson mean field theory-^ ' supplemented by phe- 
nomenological Fermi liquid parameters to account for the 
low energy quasiparticle properties''. In particular the 
superfluid stiffness and the critical temperature was cal- 
culated within two different models of the inhomogeneity. 
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The first model described doping variations on meso- 
scopic scales, comparable to or larger than the super- 
conducting coherence length. Technically we computed 
the transverse electromagnetic response tensor of a model 
system with metallic stripes embedded in an underdoped 
superconducting bulk. This was done by solving the 
appropriate Bogoliubov-dc Gennes equations within the 
rcnormalizcd mean field theory. We then averaged over 
the different stripe orientations to obtain an isotropic su- 
perfluid stiffness. 

In the second model we considered microscopic impu- 
rities that carried an excess doping charge. The tempera- 
ture dependent stiffness in this case was calculated using 
a perturbativc expansion expansion assuming both dilute 
and weak impurities (Born approximation). 

In both models, the regions of higher doping add to 
the total carrier density and hence increase the super- 
fluid stiffness at zero temperature. On the other hand 
the impurity regions give rise to low energy states that 
lead to a faster reduction of the superfluid stiffness with 
temperature. Nevertheless we found that for a range of 
doping levels in the underdoped regime a higher Tc than 
a uniform superconductor of the same doping can be at- 
tained. Moreover, in the case of microscopic impurities 
it is even possible to attain a higher critical temperature 
than the maximal Tc obtained in the pure system, that 
is, higher than Tc at optimal doping. 

The last result can help to understand the enhance- 
ment of Tc seen at the interface between an underdoped 
and a highly overdoped LSCO film.^' We have previ- 
ously noted that such an increase in Tc due to coupling 
between two homogeneous layers with different doping 
requires unrealistically strong coupling between the two 
CuO planes.^'' However if, due to the structure of the in- 
terface, overdoped and underdoped layers interpenetrate 
each other, then the proximity coupling can be induced 
by the strong in-plane hopping and the situation becomes 
equivalent to the one considered here. 

Finally we remark that our main result for the case 
of microscopic impurities was obtained within a pertur- 
bativc expansion in the impurity strength. It would be 
interesting to compare this to a numerical solution that 
takes into account scattering, at least from individual im- 
purities, exactly. If indeed excess dopants concentrated 
at random locations can lead to increase of the maxi- 
mal Tc, this opens up intriguing possibilities for further 
enhancement of Tc- For example through design of an op- 
timal ordered arrangement of the highly doped regions. 



Appendix A: The self energy in Born approximation 

We write down the Dyson equation for the disorder 
averaged Green's function 



(Al) 



From the Dyson equation we obtain the disorder averaged 
self energy, up to second order in the disorder potential 

Z^kk' = J^kk'O'S + Vkk'CTi, 

i;(k,iLj„) = nj[Z^kk + y^^kk'G'k\cu„^k'k]- (A2) 



Using (A2) we can now calculate the Nambu components 



of i;(k, Un) = Eq=o ^a^o- 



U 

So(k, w„) = -iUnUi 



with E]^ = + A|]i/2. In the hmit of small w„ this 
becomes 



So(k, w„) 



6a 
Ao 



2 



toff 



(^k + 2fiy 



SoiuJn) 

(A3) 



where Si = lim^„_^o Ek ^k/('^n + K) - V^^a/tcff/n. 
In (A3) we used the fact that in the limit of a;„— ?>0, Si = 
Ek + ED^l- Ek ^IK^I + El)- The sum 

is logarithmically divergent in the w„— )-0 limit. 



So 



E 



1 



Lot, 



log 



(A4) 



Where C is an upper cutoff for the momentum sum. To 
solve for the zero frequency limit of the self energy we 
follow Ref. -34 and assume a self consistent solution of 
the form I]o(k, w— >-0)— >■ — iFk. For the self consistent 
solution we perform the analytic continuation za;„— >w+t(5 
and replace w by its rcnormalizcd value a) = cj — S— s-iFk. 
This gives the following equation for Fk, 



1 



m 4:TTVfVA 



log r2 



C2 

2 
k 



Ao 



^1 



tcS 



a -Si) 
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where we denote 



(Ck + 

In the limit k— >-k„ode we approximate C7^ ~ fi^Sf/t'^g. 
Finally, using the fact that | {Sa/Aq) Si + {dt/tos) {1 - 
•Si ) I < 1 <C l/rii we obtain the low frequency long wave- 
length limit of the decay rate F, 



(A5) 
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Where C ~ Aq is a high energy cutoff. We shall now 
estimate the exponent and show that F is negligible. We 
plug in the doping dependent values ^ cx pit,tcs_/5t = 
Vf I[2y/2{g{p2) - g[pi))t],VAlv} ^ - 'ipi) and Ui < 
0.25. We obtain T < Aoe-i°° ~ 0. 

Using the fact that Sq— ''O, we can calculate the other 
two components of the self energy, 



k' 



E3(k,a;„) = rii?7kk - 



^'(^kk' - ^kk') - 2Ak^[/kk'lik' 

w2 + El,, 



The object Gkk' is defined by 

Gkk' = G^^kk' + G^Pk'-k^kk'Gk' + 

+ E^kPp-k^kpGpPk'-p^pk'Gk' 



As usual, the average over realizations amounts to inte- 
grating over all possible impurity positions. In all the 
summations, the only dependence on impurity positions 
appears in factors of pk- The disorder averaging gives 



Pk 



-ikr 



-ikrip — ik'r 



= ?^^^5k 

?^^'5k+k' 



(B3) 



We perform the momentum summations in the a;„— >0, 
and express the results in terms of Si as in the case of 
So. This gives Eqns. (24). 



Appendix B: Diamagnetic response 

The diamagnetic response stems from the second order 
term in the vectors potential. 



(Bl) 



-i^A,(q)i^.(-q-q')A,(q') 



q.q' 



where ^^.(r) = - J2r,atxi'^)(.'i+x,a(^r,cr + h.c. Its Fourier 
transform to momentum space is then with 



(-q-q')^Ee'^''^'''''"^-W 

r 

= J2 t{r)e^'i+'i'+^-^'^'' (^e'''^ 



(B2) 



r,k,k .cr 



Performing the disorder average leads to a diagrammatic 
expansion with three contributions, 



(i^,(-q-q')) = D0 + D1 + D2 

where Dq is the diamagnetic contribution including only 
self-energy corrections to the Green's function, described 
by diagram 5(a). An explicit calculation of this diagram 
gives 



k.k' 



e-'^^ tr {GktTajcJ, 



kk' 



Note that in our notations the trace includes the Mat- 
subara summation and Nambu space tracing. The con- 
tribution Di + D2 of diagrams 5(b) and (c) is due to 
modification of the hopping on the impurity sites and is 
therefore proportional to (5j , 

Di+D2= 5t {e''-" + e-'"'^) tr {pq+q'+k-k'Gkk'^a} 

k,k' 



We perform the sums and and keep terms up to first 
order in n,; and second order in the disorder strength dt 
and (5a- This gives 



Dq = Sq+q'2gty^cosk^ n(|^k,Ak) 

k 

Di (5q+q'4rii(5j ^cosfc^ ri(^k, Ak) 



(B4) 



D2 



kk' 



tr {GS:Zikk'Gl^,a3} 



where n(^. A) = [1 - ^/(VC^ + A2)]. Note that ^ and 
A appear in their rcnormalized values in Dq and are un- 
rcnormalizcd in Di . 

Since D2 is non vanishing at T = 0, it does not de- 
pend necessarily on low energy quasiparticles. Indeed, 
it includes a sum over all occupied states. We evalu- 
ate it numerically and express the result as a function of 
v = Ao/2toff = VA/vf in the limit of half filling. Any 
deviation from half filling introduces a small value of the 
chemical potential /i which we neglect in this calculation. 
The result has the general form 



Do 



^Pi(.) + Ap,(.) 

^0 Icff 



(B5) 



and turns out to give a negligible numerical contribution 
in the relevant regime of parameters. 



Appendix C: Paramagnetic response 

the paramagnetic current-current correlator for a given 
disorder realization is 

n,,(g,<j',c.) = - / dre'-"-(j,(g,T)j,(-g',0)). (CI) 

Naturally, after disorder averaging all contributions are 
proportional to ^qq'. The current operator is modified 
by the disordered hopping, and has the form j(q) = 
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j('^^(q) + 5j(q). The uniform part of the current operator 
in the x direction is 

= >/kqcj(.^Ck-q,CT- 

The disorder contribution to the current operator is given 
by 

6j{<D = z6i J2 Pq-k+k'(e"^- - e-'<)4^ckv. (C2) 

As a resuh, the current-current correlator has the form 
0.(q,r)j.(-q',0)) ^ (ji°)(q, r)j ^ (-q', 0)) 



gt 



^Pq_kO-i°Hk,r)4o)(-q',0)) 



, Pq-kPk'-q 



,0-(°)(k,r)4°)(-k',0)). (C3) 



fcfc' 



The disorder averaging of n^^iQi q\ ^) amounts to av- 
eraging the correlator (C3) over realizations. 

The first line of (C3) corresponds to diagrams Ho 4- 
IIi in Fig. 5, which incorporate the effects of self energy 
rcnormalization and standard vertex corrections, 



Jo 

The disorder averaged correlator takes the form 



K<i' ^kq Jk',-qtr {Gk'+q,k(-T)G'k-q,k'(T)}. (C4) 



kk' 



The disorder averaged Green's function product has a 
vertex correction part Hi which vanishes, as we will show 
below. As a result we are left with a simple product of 
disorder averaged Green's function. 



Tin 



(C5) 



J^lJkqPtr {Gk' + q,k(-T) Gk-q,k'(r)} 



= - dre'^-^^qq, V|Jfe,|2tr {Gk(-T)Gk-q(T)}. 

•^0 k 

To calculate Hq wc notice that the disorder averaged 
Green's function Gk,tj„ differ from the bare one G^ by 
renormalizcd values of Ak and ^k, as specified in (24), 



leading to renormalizcd spectrum parameters Vf and 
according to (25). This gives the result shown in (27). 

In order to sec that the standard vertex correction IIi 
vanishes, we write it explicitly as 



Ui=ni dre'""^y^ JkqJk'.-q 
Jo ,,,,, 



kk' 

Xtr {Gk'+qWk'+q,kGkGk_qZ-/k-q,k'Gk,} (C6) 

In the limit of q— >'0, the sum over momenta becomes 
4nit^ j^, sin fc^ sin fc^J^(k, k'), where is symmetric 
with respect to k and k'. Thus, under the summation 
over k or k', this contribution vanishes. 

The second line of (C3) corresponds to diagrams 112 + 
lis and the third line to 114. These contributions do not 
appear in the case of standard on-site disorder because 
they stem from direct rcnormalization of the current op- 
erator Jx(q) by an amount proportional to S^. The first 
part of this contribution is 

n2 + Ha = - / dre'^-^-^ ^ JkpJk',-,' 







kk' p 



X Pq_ptr {Gk'+q',k(-T)Gk-p,k'(T)}- 

When we insert (B3) and perform the disorder average 
we obtain a contribution of 0{U), 

n2(q,0) = -4n,^i^|Jkqptr{G^_„„G°_q,„J 

" ' k,n 

and a contribution of 0{U^), 



n3(q,0) 



^'^'T;; « -^kpJk-q -q 



9tf3 



kp.n 



tl" {Gk,LJ„Gk-p,LJ„^k-p,k-qGk_q^^„}- 



Taking the trace and performing the summations at the 
limit q— i'O we obtain the results in Eqn. (30). 

Finally, the third line of (C3), corresponding to sec- 
ond order corrections of the current-current correlation, 
yields the sum 



n4(q,0) = -(^) 1 ^ JkpJk',-p' 



9tJ (3 



pp'kk' ,n 



'+p',k,w„ Gk— p,k',ti;„ } 



5f\^ 1 V- 



pp kk ,n 



I'/kphtr {Gk^^Gk_p 



After summation this gives Eqn. (29). 
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